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We show that information about scattering data of a quantum field theory can be obtained 
from studying the system at finite density and low temperatures. In particular we consider models 
formulated on the lattice which can be exactly dualized to theories of conserved charge fluxes on 
lattice links. Apart from eliminating the complex action problem at nonzero chemical potential p, 
these dualizations allow for a particle world line interpretation of the dual fluxes from which one 
can extract data about the 2-particle wave function. As an example we perform dual Monte Carlo 
simulations of the 2-dimensional 0(3) model at nonzero p and finite volume, whose non-perturbative 
spectrum consists of a massive triplet of particles. At nonzero p particles are induced in the system, 
which at sufficiently low temperature give rise to sectors of fixed particle number. We show that the 
scattering phase shifts can be obtained either from the critical chemical potential values separating 
the sectors or directly from the wave function in the 2-particle sector. We find that both methods 
give excellent agreement with the exact result. We discuss the applicability and generality of the 
new approaches. 


In many quantum field theories the low energy exci¬ 
tations are very different from the field content of the 
elementary Hamiltonian or Lagrangian. Indeed excita¬ 
tions, rather than fundamental degrees of freedom, are re¬ 
sponsible for many important phenomena, ranging from 
superconductivity in condensed matter to confinement 
and dynamical mass generation in Quantum Chromody¬ 
namics (QCD). These excitations are often referred to 
as particles (or quasi-particles in condensed matter liter¬ 
ature) with particular masses, charges and interactions, 
although their internal structure may be rather compli¬ 
cated. For a satisfactory theoretical description non- 
perturbative methods are needed, often complemented 
with numerical simulations. 

In the approach discussed here we consider lattice field 
theories with a chemical potential p. In recent years there 
was quite some progress with finding so-called dual repre¬ 
sentations (see mm for reviews), which provide an exact 
mapping to new (dual) degrees of freedom, which are 
fluxes for matter and surfaces for gauge fields. Initially 
dualization was motivated to overcome the complex ac¬ 
tion problem (at nonzero p the action S is complex and 
the Boltzmann factor e~^ cannot be used as a probabil¬ 
ity in a Monte Carlo simulation). Here we show that very 
precise information about scattering phases and interac¬ 
tions can be obtained based on two simple observations: 

1. At finite volume and zero temperature the one, two, 
three, etc. charge sectors are separated by finite 
energy steps, and hence these charges will appear 
in the system at certain chemical potential thresh¬ 
olds pi, p 2 , P 3 , ■■■ the difference of which contains 
information about the interaction energy of parti¬ 
cles carrying charges. 


2. The dual variables allow a particle world line in¬ 
terpretation from which the multi-particle ground 
state wave function can be obtained. The latter 
in principle contains all the information about the 
interactions. 

It is these two observations which are the cornerstone of 
the two methods we develop in this work: the charge con¬ 
densation method and the dual wave function method, 
respectively. 

We develop the ideas using the 2-dimensional 0(3) 
model, and discuss their generality in the end. In ad¬ 
dition to being an important condensed matter system, 
the 0(3) model is a widely used toy theory for QCD be¬ 
cause of common features such as asymptotic freedom, a 
dynamically generated mass gap, and topological excita¬ 
tions. Moreover, there exist exact results [3] for the scat¬ 
tering phase shifts, i.e., the physical observables which 
we aim at in our new approach. Suitable dual represen¬ 
tations at nonzero p are known mm and basics of the 
finite density behavior were discussed for the related case 
of the 0(2) model in [B]. 

The Lagrangian of a quantum field theory typically 
shows some continuous symmetries resulting in corre¬ 
sponding conserved Noether charges. These symmetries 
and charges play a two-fold role in our approach: When 
transforming the theory to the dual representation the 
symmetries give rise to constraints of the dual variables 
making flux conservation explicit. The resulting “world 
lines of flux” transport the corresponding charge in space 
and time, and the particle number turns into a topolog¬ 
ical invariant: the winding number of the corresponding 
flux around the compact Euclidean time. 

The second important aspect of the symmetries and 
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the corresponding charges is that a chemical potential 
can be coupled to the associated charge, which then can 
be used to populate the corresponding charge sector. In 
the dual representation the chemical potentials introduce 
an additional weight for the winding number of the cor¬ 
responding flux lines. 

At sufficiently low temperature T and finite volume 
one can control the particle number with the chemical 
potential and systematically probe the system in the dif¬ 
ferent charge sectors, which are separated by the afore¬ 
mentioned critical values ^i,i = 1,2,3... of the chemical 
potential. The critical values are clearly marked by 
steps when plotting the expectation value of the charge 
Q at low T as a function of /r (see [6] and Fig. [^below). 
In the presence of a mass gap m, the (‘Silver blaze’) range 
from /r = 0 to /i = p,i corresponds to the charge-0 sector, 
no particle is present, and for sufficiently large spatial 
volume one has ni = m. The interval /i S (/ii,/i 2 ) de¬ 
limits the charge-1 sector, where ^2 corresponds to the 
energy of the lowest charge-2 state, i.e., containing two 
particle masses plus their interaction energy [T3j. At fi¬ 
nite volume, it is exactly this 2-particle energy that can 
be related to scattering data using the Liischer formula 
[ZHS], and analyzing the 2-particle condensation thresh¬ 
old /r 2 as a function of the volume constitutes our first 
approach to extract scattering data from nonzero den¬ 
sity and finite volume. We refer to this approach as the 
“charge condensation method”. 

Our second approach - which we refer to as the “dual 
wave function method” - is based on a direct analysis of 
the dual fluxes in the 2-particle sector, i.e., the interval 
/i G {gL 2 ,gLs), where /is marks the onset of the charge- 
3 sector. In this interval we analyze the two winding 
flux loops that characterize the charge-2 sector and de¬ 
termine the distribution of their spatial distance which 
we relate to the 2-particle wave function, from which we 
again compute scattering data. 

Both methods are presented in detail for the 2-D 0(3) 
model. We begin with discussing the dual representation 
as derived in [?] where the details of the dualization and 
also the conventional form of the model are presented. 
We consider the 0(3) model with a single chemical po¬ 
tential coupled to one of the conserved charges. The 
conventional degrees of freedom are 0(3) rotors and the 
chemical potential is coupled to the 3-component of the 
corresponding angular momentum. In the dual form the 
partition sum Z is exactly rewritten into a sum over con¬ 
figurations {m, I, k} of 3 sets of dual variables mx,u G Z, 
lx,i/,kx,v G IN on the links {x,v) of a W x Ns lattice 
(periodic boundary conditions, lattice constant a): 

Z{g) = Bj[m,l,k] ( 1 ) 

{ m, Z, fc} 

X - mx-0,S)- 

X V 

Each configuration of the variables vrix^vi comes 

with a real non-negative weight I, k] which depends 



FIG. 1: The charge Q as a. function of g (in units of the mass 
m). We compare different temperatures T at fixed spatial 
extent Lm = 4.4 {Nt = 100, 200,400,1000 at Ns = 20 using 
J= 1.3). 


on the coupling J. Bj[m,l.,k] is given explicitly in [3], 
but irrelevant for the discussion here. A second weight 
factor originates from the chemical potential 

/i which couples to the temporal component of the cur¬ 
rent mx,v (in our notation ^ = 0,1 and v = 0 denotes 
Euclidean time). Obviously all weights are real and pos¬ 
itive such that the complex action problem is solved in 
the dual formulation for arbitrary /r. 

The dual variables mx,u, which correspond to one of 
the 0(2) subgroups of 0(3), obey constraints at each 
site X. The product of Kronecker deltas i5(..) enforces 
'S/fhx = — 'mx-i/,A = 0 Vx, which is a discrete 

version of a vanishing divergence condition. Thus the 
flux mx^u is conserved and the admissible configurations 
of TOj,-flux are closed loops. 

Since the TOi/-flux must form loops, the term in Q that 
multiplies /r can be written as a^x’^x.o = aNtw[mi,], 
where w[7nu\ is the total winding number of ruiz-flux 
around the compact time direction with extent aNt = 
I/T (we use natural units with ks = h = c = 1). Thus 
we identify the winding number w[m^] as the particle 
number in the dual formulation. 

The constrained m^^-fluxes can be updated with a gen¬ 
eralization of the worm algorithm nni, while for the other 
fluxes local Monte Carlo updates are sufficient. Most 
of the Monte Carlo results presented in this letter were 
computed at fixed coupling J = 1.3 (with some scaling 
checks performed also at J = 1.4 and J = 1.5, i.e., closer 
to the continuum). The temperature T = 1/aNt was 
varied by changing the temporal extent W of the lattice, 
the spatial size L = aNg by changing Ng. Dimension¬ 
ful quantities are expressed in units of the mass m of 
the lowest excitation, which was determined from prop¬ 
agators in the conventional representation, and at, e.g., 
J = 1.3 is am = 0.222. 

To substantiate the above discussion of sectors with 
fixed Q, we show in Fig. our J = 1.3 results for 
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FIG. 2: Left: The susceptibility x at fixed temperature T/m = 0.0045 as a function of chemical potential for different spatial 
extensions L {Nt = 1000 with Na = 20,30,40 and 50). For each L we show the first two peaks corresponding to the critical 
chemical potentials ni and 112 ■ Right: The critical chemical potentials at the same temperature now as a function of L. The 
solid curves are results of fits (exponential for fii and power-law for ^ 2 )- 


Q = {w[mi,]) as function of /r for several low temper¬ 
atures, at fixed mL = 4.4. Decreasing the tempera¬ 
ture we indeed find the expected formation of plateaus 
in Fig.[^ which correspond to the sectors of fixed charge, 
and we can read off the critical values In a practi¬ 
cal calculation one actually determines the fj,i from the 
peaks of the corresponding particle number susceptibility 
X = ((w[mjy] — QY)IL. These susceptibilities are shown 
as a function of /r in the left panel of Fig.[^ and we com¬ 
pare results for different spatial extents L in units of to. 
The susceptibilities show pronounced peaks which we can 
use to determine the values We remark that the po¬ 
sition /ii of the first peak is independent of L, while the 
second peak ^2 shifts to smaller values when increasing 
L (see the discussion below). 

We now make quantitative the above arguments con¬ 
necting the critical chemical potentials /ii and \i 2 with 
the mass to of the lightest particle and the energy of the 
2-particle states. We write the grand-canonical partition 
sum .^(/i) with a grand potential D(/i) in the form (iJ 
and Q denote Hamiltonian and charge operator) 

Z{^i) = tr = g-0(M)/T^ (-2) 

The small-T limit is governed by the minimal exponents, 
i.e., in each sector with charge Q the corresponding min¬ 
imal energy dominates and the grand potential in 

the different sectors is 

for/r e [0,/ii) , 

I - fi = m- for /r e (/ii, /ia) , 

I - 2fi = W -2^J. for /r e (/i 2 , /is) , 

( 3 ) 

where we introduced the minimal 2-particle energy W. 
Note that in the charge-1 sector we assumed fj,i = to, 
i.e., we neglected possible finite volume corrections (see 
below for a cross check). 


The 2-particle energy W can be calculated from 
and ^,2'- The first two transitions between neighboring 
sectors occur when 0 = to — ^1 and to — /i 2 = FF — 2 /i 2 , 
respectively, and thus we find 

IF = TO -I- /i2 = /il + Ai2 . (4) 

So far we have not discussed the role of nonzero spa¬ 
tial momenta of the states, which for our 2-dimensional 
model are given by 27rn/L = 2TTn/aNs, n = 0,1... Ns—1. 
We stress again that at very low T only states with van¬ 
ishing total momentum contribute to the partition sum. 
Thus in the charge-1 sector we have one particle at rest, 
in the charge-2 sector two particles with opposite spa¬ 
tial momentum et cetera. We can combine this with 
the discreteness of the momenta to understand the in¬ 
dependence of the fii, which is documented in the right 
panel of Fig. where we plot the values of fii and ^2 as 
a function of L. 

For Hi we expect a dependence on L only when L be¬ 
comes smaller than the Compton wave length of the light¬ 
est excitation m and self interactions around periodic 
space alter the mass. This is indeed what we observe in 
the right panel of Fig. For H2 a-nonvanishing depen¬ 
dence on L is expected throughout, since the box size 
L controls the allowed relative momentum. Also such a 
non-trivial dependence is obvious from Fig.[^ 

For a short ranged potential one can write the two par¬ 
ticle energy W as twice the energy of an asymptotically 
free particle, 

IF = 2\/to 2 -b , (5) 

where k denotes the relative momentum which is shifted 
from the values 21:11 jL of the free case. In a finite vol¬ 
ume, only certain quantized values of k can account for 
the scattering phase shift S{k) of the interaction and the 
periodicity, as expressed by the Liischer formula [7] 

^ 2 iS{k) _ g-ikL 


( 6 ) 
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FIG. 3: Results from the charge condensation method (blue 
symbols, labelled as CC) and the dual wave function method 
(red, DWF) using different J and Ns- We compare our results 
to the analytic solution [3] (full curves) and the analysis based 
on 2-particle spectroscopy [9] (blue, LW). In the top panel we 
show the momenta fc as a function of L and in the bottom 
panel the phase shift S{k) versus k {k and L in units of m). 


Varying L allows one to scan a whole range of momenta 
k. The temperature must be low enough for pronounced 
plateaus to form (see Fig. [^, which gives rise to the fol¬ 
lowing two conditions: T m and T/m ^ l/(Lm)^. 
Our numerical results for this extraction are given in 
Fig. [31 In the top panel we show our data for fc as a func¬ 
tion of L and compare to the exact result [5] for “isospin 
2”, which is the relevant case for our choice of the chem¬ 
ical potential which excites the 3-component of the 0(3) 
angular momentum. We also include results from the nu¬ 
merical spectroscopy calculation in the 2-particle channel 
[9] and find excellent agreement of our results with the 
analytical and numerical reference data. In the bottom 
panel we give the results for the phase shift 6{k) as a 
function of fc, and again find excellent agreement with 
the reference data. 

Our second approach for the determination of scatter¬ 
ing data, the “dual wave function method”, is based on a 
direct analysis of the flux variables which carry the 
charge: In the charge-2 sector we identify the two flux 
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FIG. 4: The relative wave function iIj{Ax) as determined 
from analyzing the flux lines in the charge-2 sector. We 
show results for J = 1.3 at T/m = 0.0045 (W = 1000) and 
different spatial extents, Lm = 4.4,6.6,8.8,11,13.2 {Ns = 
20,30,40,50,60). The wave functions are fit according to 
Eq. 0 and in the legend we give the corresponding fit re¬ 
sults for fc. 

lines which wind around the compact time and interpret 
their temporal flux segments = 1 as the 

spatial position a;i of the charge at time xq. Thus for 
a given time Xq we obtain two positions x^^'^ and 
identify Ax = — a;j^^)| as the distance of the two 

charges at that time. Sampling over time and many con¬ 
figurations we obtain the probability distribution of the 
distance Ax, and the square root of this distribution can 
be identified as the relative wave function 'tjj{Ax) of the 
two charges. 

In Fig. [^ we show 'ip{Ax) for different spatial extents 
L. Except for very small Ax the wave functions are very 
well described by shifted cosines, 

'ip{Ax) oc cos{k{Ax — L/2)) . (7) 

This finding confirms the applicability of Eq. (j^ : Outside 
the interaction range the wave functions for the relative 
motion of the two particles are standing waves with mo¬ 
menta fc which are related to the scattering phase shifts 
via §. Thus we can fit the data for ^(Ax) with the 
cosines 0 and obtain the momenta fc from that fit. The 
fit results are also shown in Fig. [^ and in the legend we 
give the corresponding momenta fc. From these one again 
obtains the phase shift 6{k) via §. The results from the 
dual wave function method are included in Fig. [^ and 
agree very well with the charge condensation method and 
the reference data. 

Let us summarize the two approaches for obtaining 
scattering data presented here and comment on their 
generality: Both methods are based on simulations at 
nonzero chemical potential /i at very low temperatures. 
The charge condensation method relates the critical 
chemical potentials /Xi and /i 2 to the 2-particle energy, 
which in turn can be related to the scattering phase shift 
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via the Liischer formula. Technically this method as¬ 
sumes that nonzero chemical potential simulations are 
feasible, which for example is the case for the isospin po¬ 
tential in QCD. The basic concept of the charge conden¬ 
sation method is generalizable also to higher dimensions, 
in particular the relation of the condensation thresholds 
/ii, fi 2 to the 2-particle energy. In higher dimensions a 
full partial wave analysis would be necessary for complete 
information about scattering, but at least the scattering 
length can be extracted from the 2-particle energy [Hi. 
For the second approach, the dual wave function method, 
the chemical potential is adjusted such that the system is 
in the charge-2 sector. This method is rather general in 
arbitrary dimensions. The relative 2-particle wave func¬ 
tion can be determined from the fluxes in the charge-2 
sector and in principle this gives access to the complete 
scattering information. 

We remark that both methods can be straightfor¬ 
wardly generalized to three and more particles. A second 
remark concerns the possibility to couple different chem¬ 
ical potentials to different conserved charges. Using such 
a setting one can also populate 2-particle sectors with 
two different particle species and study their scattering 
properties. While the charge condensation method is in¬ 
dependent of the representation of the system, we em¬ 


phasize that the dual wave function method requires a 
suitable dualization and makes explicit use of the word¬ 
line interpretation. The latter method opens the possi¬ 
bility to study particles and antiparticles, i.e., charges of 
opposite sign. 

This exploratory paper is a first step towards a full 
understanding of the relation between low temperature 
multi-particle sectors and scattering data. We believe 
that this is an interesting connection to explore, and un¬ 
derstanding the physics involved clearly goes beyond the 
development of a new method for determining scattering 
data in lattice simulations. 
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